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SUMMARY 

^ A brief  description  is  given  of  a ray  tracing  nodel  developed 
for  Baking  propagation  loss  predictions  for  use  in  underwater 
acoustic  research.  A range  independent  sound  speed  profile 
and  ocean  depth  are  assuned.  The  various  options  available  to 
the  user  are  discussed  together  with  some  results,  which 
highlight  the  iaportance  of  having  available  good  bottom  loss 
versus  grasing  angle  data.  A FORTRAR  computer  programme  is 
currently  available  for  use  on  the  IBI  370/168  at  Defence 
Research  Centre  Salisbury. 
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1 . IHTBODUCTIOH 

This  paper  describes  the  operation  of,  and  sobs  results  fron  a computer 
programme  used  in  propagation  loss  calculations  for  underwater  acoustic  research. 
Ray  theoretical  expressions  are  employed  in  the  calculations  and  an  ocean  with 
a flat  bottom  and  a range  Independent  sound  speed  profile  are  assumed.  The 
main  part  of  the  programme  was  written  by  Dr.  H.  Sucker  of  the  Haval  Ocean 
Systems  Center  (formerly  the  Haval  Undersea  Center),  San  Diem.  The  programme 
was  adapted  for  use  on  the  IBM  7090  and  later  on  the  IBM  370/168  at  WSRL  and 
has  been  in  use  for  several  years.  During  this  time  a number  of  modifications 
and  additions  have  been  made,  which  are  described  in  this  paper. 

The  programme  was  written  to  provide  fast  and  accurate  propogation  loss 
estimates  from  a source  to  a number  or  receivers.  It  employs  a ray  sweep  out 
method,  which  is  described  in  Reference  1 . This  method  identifies  all  signifi- 
cant  ray  paths  from  a source  to  a receiver  during  a single  sweep  out  of  the 
source  beam  patterns  and  at  the  same  time  calculates  the  total  ray  intensity  at 
the  receiver. 

Ho  attempt  is  made  to  describe  the  details  of  the  programming  techniques  used. 
However,  the  input  data  required  to  use  the  program  on  the  DRC  1^  370/l68  are 
Included  in  Appendix  I.  The  names  of  the  more  Important  input  variables  are 
included  as  bracketed  capitals,  when  the  various  options  are  being  discussed  in 
the  text. 


2.  APFROXINATIHG  THE  SOUHD  SPEED  PROFILE 

The  sound  speed  profile  used  as  input  to  this  and  many  other  propagation  loss 
programmes  \isually  consists  of  set  of  sound  speed  values  at  a number  of  specified 
depths.  These  values  say  have  been  read  off  a continuous  trace  of  a sound  speed 
profile  or  be  the  results  of  measurements  at  discrete  depths  (e.g. , resulting 
from  Hansen  cast  measurements  of  temperature  and  salinity).  If  the  depth/sound 
speed  values  have  been  obtained  from  a continuous  trace  then  the  values  are 
selected  so  that  the  important  characteristics  of  the  sound  SMed  profile  are 
retained  (e.g.  depth  of  mixed  layer,  depth  of  SOFAR  axis  etc.).  It  is  important 
to  keep  tte  number  of  depth/sound  speed  values  fairly  small,  since  the  ccmputing 
time  increases  with  the  number  of  values  employed;  ten  values  are  usually 
adequate  unless  the  profile  is  imduly  complicated. 

Given  a set  of  depth/sound  speed  values  the  programme  offers  three  options  for 
fitting  a profile.  A brief  description  will  be  given  of  these  options. 

2.1  Linear  interpolation 

The  first  is  a linear  interpolation  (ITTFKO  » O)  between  the  input  deptl^ 
sound  speed  values.  The  profile  is  linear  in  the  variable  l/0*,  where  C 
is  the  local  sound  speed  and  ■ is  the  depth; 

l/C*  ■ a ♦ ba,  ....  (l ) 

where  a and  b are  oonstants.  Use  of  this  expression  la  common  praotiee  in 
aeouetios,  sinoe  it  allows  exact  solutions  to  be  found  for  the  range 
travelled  by  a ray  over  a given  depth  interval  (e.g.  Appendix  II  and 
referenee  2).  Seleotion  of  this  option  results  in  a profile  made  up  of 
sei^nta  satisfying  equation  (l)  and  paaaing  through  the  given  depth/aound 
apeed  values. 
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2.2  natural  apline 

In  ordar  to  eliainate  apurioua  Intenaity  values , which  arise  when  the 
gradient  of  the  sound  speed  profile  is  diseontinuous(raf.3)»  a natural  spljne 
can  be  fitted  through  the  soiuid  speed/depth  values  using  the  second  of  the 
profile  fitting  options.  The  natural  spline  is  a cubic  polynomial; 

l/C*  ■ a bs  + ca*  + ds®,  ....  (2) 


which  has  continuous  first  and  second  derivatives,  passes  throtjgh  the  input  J 

data  values  and  has  zero  valued  second  derivatives  at  the  boundaries.  The  | 

coefficients  (a,  b,  c and  d)  are  determined  by  an  iterative  numerical  pro- 
cedure for  each  segment  of  the  profile.  Care  should  always  be  exercised  in 
the  selection  of  the  input  data  values  so  that  the  fitted  curve  is  a reason- 
able replica  of  the  original  profile.  This  can  be  a problem  where  there  are 
sudden  changes  in  the  sound  speed/depth  profile,  e.g.  at  the  bottom  of  the 
mixed  layer  or  if  step  changes  occxir  in  the  themocline.  The  curve  fitting 
procedure  sometimes  results  in  overshoot  and  an  unacceptable  orofile.  Use 
of  alternative  methods  for  calculating  the  spline  (e.g.  ref, 4)  may  improve 
the  fit;  these  have  not  been  investigated.  ^ 

2.3  Natural  spline  with  least-squares  fit  . 

The  third  option  available  requires  two  sets  of  sound  speed/depth  data.  j j, 

A spline  is  fitted  to  the  second  set  and  the  mean-square  difference  between  ' 

values  of  l/C*  determined  from  the  spline  and  those  from  the  first  data  set 
is  computed.  The  l/c^  values  in  the  secoxtd  set  are  then  modified,  a new 
spline  fitted  and  the  mean-square  difference  evaluated.  This  process 

continxies  until  the  mean-square  difference  is  less  than  a specified  value.  | 

This  option  has  not  been  used  in  aoy  of  the  studies  using  the  programme,  but  !. 

may  be  useful  if  severed  sets  of  data  are  available  over  an  area  of  interest  j 

and  a best  fit  needs  to  be  found.  | 


3.  SURTACE  INAOIHG  EFFECTS 

It  is  well  known(zef.5)  that  when  an  acoustic  source  is  near  the  surface 
(i.e.  within  several  wavelengths)  the  propagation  loss  to  a distant  receiver  is 
increased  over  that  for  a deeper  source.  This  arises  because  of  interference 
between  the  surface  reflected  and  direct  ray  paths  and  is  sometimes  referred  to 
as  surface  decoupling.  The  effect  is  of  special  interest  in  calculating  the  j | 

radiated  noise  from  surface  ships,  since  they  radiate  principally  at  low  fre-  : 

quencies  and  have  an  acoustic  centre  with  a depth  of  only  a fraction  of  a wave-  f 

length  at  these  frequencies.  | 

The  programme  does  not  include  the  facility  of  coherently  sunning  the  oontrlb-  i 

utions  from  different  ray  paths.  The  Intensity  at  a receiver  is  conputed  from 
the  spacing  of  adjacent  rays,  lAioh  bracket  a receiver  (Appendix  II ) and  the  I 

total  intensity  found  by  sunning  the  Intensities  due  to  all  ray  paira.  In  . 

order  to  include  the  near-surface  effeota  a near-surface  source  is  replaced  by  a ‘ ! 

single  source  at  the  surface  with  a dipole  radiation  pattern  (lIVBP  •3); 

P(t)  - 4 sia(2v  d/c.f  sini  ) (3)  i 

( 

where  9 is  the  ray  elevation  angle  naasurad  fron  the  horisontal,  f is  the 

frequency,  d is  the  source  depth  and  C is  the  sound  speed  at  the  surface.  1 

This  is  only  an  approzinatiom  to  the  real  situation  but  should  be  quite  good 


i 

v 
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when  the  near-surface  water  has  a constant  sound  speed  at  least  to  the  source 
depth. 

A near-surface  receiver  is  treated  in  a slifhtly  different  fashion  and  this  is 
discussed  in  the  next  section. 


4.  RAY  ARRIVAL  SISUCTllRS  AID  VERTICAL  BBAM  PATTBRHS 

As  well  as  providinj^  the  total  intensity  at  a receiver,  using  the  technique 
outlined  in  ApMndiz  II,  the  progran  can  provide  the  vertical  arrival  structure 
at  a receiver  (IJCFR  Negative).  Up  to  100  angle  increnents  ccm  be  eaployed  over 
the  vertical  angle  range  0*^  to  180”  and  the  total  intensity  of  all  rays  arriving 
at  the  receiver  in  each  increoent  (added  incoherently)  can  be  detezained. 

The  above  option  is  necessary  when  it  is  required  to  coapute  the  response  of 
an  array  with  soae  directivity  in  the  vertical  plane  (iXFR  > -2).  The  vertical 
beaa  pattern  (beaa  in  dB  units)  is  provided  as  input  data.  The  nuaber  of  values 
aust  be  less  than  or  equal  to  the  nuaber  of  angle  increaents  above.  If  less,  the 
last  given  value  of  the  beaa  pattern  is  assuaed  to  bold  over  the  reaainder  of  the 
vertical  angle  range.  The  vertical  beaa  pattern  is  then  convolved  with  the 
acoustic  arrival  pattern  and  the  array  output  for  all  steered  directions  of  the 
array  detezalned. 

The  above  procedure  is  only  approxiaate,  since  it  assuaes  that  the  sound  field 
is  unlfoza  over  the  vertical  extent  of  the  array,  that  is  each  ray  is  associated 
with  a plane  wave  incident  on  the  vertical  array  and  the  total  acoustic  field  is 
the  sua  of  the  energies  of  these  acoustic  plane  waves.  The  approxiaation  will 
be  quite  good  when  the  array  is  short,  but  will  be  less  accurate  when  the  array 
extends  over  a depth  range  across  which  there  is  a significant  change  in  the 
sound  field.  This  later  case  has  to  be  treated  by  a aore  exact  method  as  is  done 
for  exaaple  in  reference  6,  where  the  sound  field  amplitude  and  phase  at  each 
receiver  in  the  array  were  computed  using  a aode  theory  model. 

It  is  possible  to  input  a vertical  beaa  oattem  for  the  source  in  addition  to 
the  two  options,  oanldirectlonal  (INVBP  > o)  and  dipole  (INVBP  * 3),  previously 
discussed.  Using  another  option  (IRVBP  > 1 ) the  vertical  pattern  is  entered  at 
1°  intervale  for  the  vertical  angle  range  0 to  160*^  and  is  assumed  to  be  syaaet- 
rical  about  0^.  The  final  option  (lIVBP  > 2)  pezaits  two  sets  of  vertical 
patterns  tc  be  entered  covering  the  range  0 to  180<>  axtd  0 to  -ISO**  at  1°  Intervals, 
where  the  angles  are  aeasuzed  from  the  horlaontal.  It  is  not  necessary  to  enter 
all  360  values;  entering  a value  equal  to  or  less  than  sero  teralnates  the  input 
list  and  automatically  sets  all  subsequent  data  points  to  the  Icust  value  read. 

Anrther  input  paraaeter  (ITILTD)  pezaits  these  vertical  beaas  to  be  steered 
above  and  below  the  horisontal.  A positive  value  steers  the  aaln  beaa  in  a 
direction  above  the  horisontal. 


3.  BOTTOM  LOSS 

A number  of  bottom  loss  options  (IABOT)  is  available  within  the  progran.  A 
perfectly  absorbing  bottom  (lABOT  - 2)  perfectly  reflecting  bottoa  (lABOT  > O) 
or  a bottoa  loss  table  given  at  1^  intervals  (lABOT  s l)  can  be  specified  in  the 
input  data.  If  the  value  of  the  input  variable  (IABOT)  lies  between  11  and  18 
a set  of  bottoa  loss  versus  grasing  angle  values  are  selected  fron  a table 
contained  in  the  subroutine,  EHLOS.  The  set  selected  depends  on  the  frequency 
and  bottom  roui^mess  (Xl  equal  to  IABOT  -10)  specified.  This  table  was  taken 
from  the  FAiCT(ref.7)  ray  tracing  nodal. 

For  frequencies  betveea  1 and  3.5  kHs  the  bottom  losses  are  those  given  by 
the  proposed  US  levy  interim  standard  bottoa  loss  curves.  liid^t  rou^iaess 
(ll  B 1 to  8)  are  available  for  this  frequency  range.  Sets  of  bottom  loss 
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values  are  also  given  for  the  frequency  Intervals  150  Hs  and  below,  130  He  to 
700  He,  700  He  to  1 kHe,  3«5  lAe  to  4*3  kHe  and  for  frequencies  greater  than 
4.3  kHe.  Five  rou^mess  options  are  available  within  each  firequency  range 

(IR  =1  to  3).  Figure  1 shows  the  bottoa  loss  versus  graelng  cmgle  curves 
for  frequencies  less  than  130  Be  for  each  roughness  value. 

Propagation  loss  is  very  sensitive  to  bottoa  loss  emd  the  selection  of  botton> 
loss  values  versus  graeing  angle  is  a particularly  difficult  problem  if  reliable 
PL  preductions  are  required. 

Propagation  loss  predictions  (figure  2)  have  been  made  using  the  prograaoe  for 
a frequency  of  100  Hz,  for  two  water  depths  (2000  a and  3000  ^ and  for  two 
bottom  roughnesses  (iR  = 1 and  4).  The  most  striking  featiire  is  the  large 
differences  between  the  curves,  even  at  quite  modest  ranges  (e.g.  22  dB  at  100 
km). 

The  results  for  the  3000  m water  depth  provide  a good  example  of  the  sensit- 
ivity of  the  PL  results  to  the  bottom  loss  data  used.  The  major  contributions 
to  the  received  energy  at  ranges  exceeding  about  23  km  arises  from  rays  having 
grazing  angles  at  the  bottom  less  than  about  13  degrees.  Over  this  angular 
range  the  BL  values  used  in  generating  the  PL  curves  differ  by  only  2 to  3 dB. 

In  order  to  make  a prediction  of  PL  accurate  to  several  dB  at  130  km  the  BL  data 
must  be  known  to  better  than  1 dB.  Data  with  this  accuracy  are  rarely  available 
in  Australian  waters,  especially  at  the  lower  frequencies. 

The  convergence  zone  peaks  in  figure  2 are  of  course  independent  of  the  BL. 

No  correction  is  made  in  the  model  for  the  intensity  of  sound  in  the  vicinity 
of  a ray  theory  caustic  as  is  done,  for  example,  in  the  FACT  model.  The  method 
used  in  the  current  program  (Appendix  II ),  which  involves  spatial  averaging, 
usually  results  in  modest  intensity  values  within  a convergence  zone,  and 
excessively  high  intensity  values  (or  low  PL)  generally  don't  occur.  The 
variations  in  the  heights  of  the  peaks  in  figure  2 are  due  to  the  finite  range 
interval  between  successive  PL  determinations. 

Results  for  the  shallower  water  depth  of  2000  m,  show  a greater  divergence  of 
PL  values  for  the  different  BL  data,  because  of  the  greater  number  of  bottoa 
bounces  required  to  reach  a given  range.  At  130  km  the  PL  results  differ  by 
19  dB  compared  to  10  dB  in  the  deeper  water.  Conseqiiently  the  accuracy  of  the 
BL  data  must  be  increased  to  between  one  quarter  one  half  dB  if  the  PL 
preductions  are  going  to  be  accurate  to  several  dB  at  this  range.  In  still 
shallower  water,  th«  problem  would  be  worse  except  that  the  ranges  of  Interest 
begin  to  be  reduced  because  of  the  overall  increase  in  the  PL. 

In  addition  to  these  problems  associated  with  BL  data  accuracy  there  is 
another  question  relating  to  the  adequacy  of  the  description  of  the  BL  in  terms 
of  simple  Raylei^  plane  wave  reflection  coefficients  at  low  frequenoies(ref.8). 
There  is  considerable  penetration  of  sound  into  the  bottoa.  Associated  with 
this  penetration  are  reflection  from  sub-bottom  interfaces,  refraction  in  the 
sub-bottom  layers,  absorption  in  the  layers,  generation  of  shear  waves  in  the 
solid  layers  and  generation  of  surface  waves  at  the  layer  interfaces,  to  mention 
only  some  of  the  bottoa  effects.  The  present  model  is  only  oapable  of  using 
simple  reflection  coefficients;  this  may  not  be  adequate  at  low  frequencies  in 
some  areas. 

A common  feature  of  low  frequency  BL  data  (using  reflection  coefficient 
interpretation  at  small  graslng  angles(ref.8)  is  a negative  loss  due  to  refract- 
ion and  possibly  refleetlons  within  the  bottM  layers.  BL  data  of  this  type 
were  obtained  in  the  north-eastern  Indian  Ooean  during  a recent  survey  conducted 
by  the  U.S.  Naval  Oceanogrm^o  Office.  The  preliminary  results  are  given  in 
reference  9. 

In  conclusion,  the  uncertainties  la  the  BL  data  moke  FL  predlotions  dlffloult 
for  a particular  area  and  reliable  BL  data  must  be  available  for  areas  of  interest. 
Often  the  differences  between  two  models  beoome  unimportant  because  of  the  basic 
uncertainties  in  the  input  data  required  by  both. 
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6.  ADDITIONAL  FUNCTIONS 

There  ere  several  addltonal  programme  functions  which  should  be  mentioned: 

(1)  A surface  reflection  loss  (SLDB)  can  be  included  but  is  independent  of 
surface  grazing  angle. 

(2)  Volume  absorption(ref .10)  is  included  if  a frequency  value  is  provided  in 
the  input  data(raBQ). 

(3)  Line-printer  plot  of  FL  versus  range  (IPRL). 


7.  CONCLUDING  RBNABKS 

The  programme  has  been  quite  useful  for  malfl ng  a variety  of  PL  and  ray  arrival 
structuire  predictions  over  the  past  few  years.  It  has  of  course  the  limitation 
of  dealing  only  with  a single  sound  speed  profile  in  a flat  bottom  ocean,  but 
should  continue  to  be  useful  for  a wide  range  of  applications.  Current  model 
developments  at  DRCS  (e.g.  ref.1l)  are  concentrating  on  models  which  can  handle 
range  dependent  ocean  depth  and  sound  speed  profiles.  These  models  have 
immediate  use  in  two  areas  of  interest,  the  first  over  the  continental  slope  and 
the  second  in  the  vicinity  of  oceanographic  fronts  and  eddies.  However,  they 
are  demanding  of  computing  time  emd  must  be  used  with  some  care,  so  that  the 
current  ray  model  should  still  find  many  uses. 
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APFKHDIX  I 

FROGRAIOIB  IHPDT  DATA 

The  data  cards  required  to  run  the  prograa  together  with  the  required 
foraata  are  given  below. 

I 


Card  lo. 

Variable 

Poraat 

Options 

\ lotea 

1 

ITTPRO 

12 

0 

Linear  segments'] 

1 

Hatural  spline  \ (g,,  section  2) 

2 

Fitted  spline  } 

HGQ 

12 

3 

Gaussian  Qiutdrature  variable  used  if 

ITTPRO  = 1 or  2 ( Section  2) 

TITLE 

10A4 

Title 

2,3  ... 

ZNL  (l) 
VNL  (Z) 

P8.1 

F8.3 

K .pe«i}  If  IWPBO  . 2 

« 

ZILD  (l) 

P8.1 

Depth 

VIL  (I) 

P8.3 

Sound  speed 

Zero  value  of  VML  or  VIL  terminates  input 
list 

*A 

GANDD 

P3.1 

Upgoing  ray  limit  at  source 

P3.1 

Ray  cmgle  increment  for  IPLRAT  s i 

GAMDD 

P3.1 

Downgoing  ray  limit  at  source 

IHVBP 

11 

0 

Vertical  beam  pattern  • ' 

Omnidirectional 

1 

Read  pattern  at  Intervals  for  0 to  180^ 

2 

Read  pattern  at  1°  intervals  for  0 to  180° 
followed  by  0 to  -180° 

3 

Dipole  source  pattern  (see  sections  3 and  4] 

ITILTD 

12 

Beam  steer  angle t positive  for  upward 
pointing  beam 

ZPLSAT 

12 

Test  function  - not  required 

3081XV 

P5.2 

Source  level  (dB  re  ref.  level) 

SLDB 

F4.2 

Surface  loss  (dB) 

lAfiOT 

12 

0 

Bottom  loss 

Perfeotly  reflecting 

1 

Requires  loss  table  at  1°  inteinrals 

2 

Perfeotly  absorbing 

1Z 

1 e X < 8 Define  bottom  rou^mess  (see 
section  5) 

lAST 

12 

lot  used 

12 

lot  used 

mm 

T10.0 

Initial  receiver  range  (metres  or  yards) 

1^1 

FI  0.0 

Range  inoremsnt  (metres  or  yards) 

IL 

F10.0 

Final  rsosiver  range  (metres  or  yards) 
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Card  No. 

Variable 

Fomat 

Options 

Notes 

IXPR 

12 

0 

1 

2 

-1 

-2 

Print  diagnostic 

Extensive  print  diagnostic 

Intensity  vs  angle  printout  at  spacing  of 
KLT  (see  Sub  TRACE)  Max  number  of  range 
pta  is  40 

Intensity  vs  angle  followed  by  convolution 
with  receiver  beam  pattern.  Max  number  of 
range  pts  40  (see  section  4) 

FRBQ 

F5.0 

Frequency  for  volume  attenuation  calculation 

EIHD 

13 

0 

1 

2 

Printout  units  - KTDS 

- KM 

- MM 

IFSL 

1 

Line  printer  plot  of  PL  vs  Range.  Max 
number  of  range  pts  is  111. 

» 

VBPR  (l) 

F4.2 

Source  beam  pattern  (dB  units)  at  1*^ 
intervals.  If  INVBP  s 1 max  number  of 
values  is  180  (symmetric).  If  INVBP  = 2 
max  number  of  values  is  360.  Negative 
value  terminates  input. 

* 

BTL  (l) 

P4.2 

Required  if  lABOT  = 1 . Bottom  loss  (dB) 
at  1*^  intervals.  Terminate  input  list 
with  negative  value  (l  must  be  greater  tban 
9). 

• 

MT  (l) 

1214 

Required  if  IXPR  = -2.  Convolved  output 
heading. 

« 

BBiM  (l) 

F4.2 

Required  if  IXPR  * -2. 

Values  specified  at  luigle  increments  of 

RLT  (Sub-routine  TRACE).  Tarmlnate  input 
list  with  negative  value. 

• 

SD 

RCI)  (l) 

F8.1 

F8.1 

• 

Source  depth 

Receiver  depths  (number  of  range  pta  x 
number  of  receivers  less  than  4000). 
Terminate  list  with  null  value. 

• 

IK»B 

15 

Von  Zaro 

0 

Read  new  input  from  *A  (i.a.  uae  axiatiag 

profile ). 

Tarmlnate  run 

■ORSt  (l)  *r«pc«MixtB  nav  eu^a) 
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APFBHDIX  II 

RAY  THEORY  RELATIONSHIFS 

Soae  of  the  rsytheory  relationehips  used  to  derive  intensity  -rnwga  values 
in  the  aodel  are  included  in  this  section. 

II. 1 So\ind  intensity 

The  intensity  range  r relative  to  the  intensity  at  unit  range 

(l^)  is  given  by  (e.g.  ref. 12,  page  57) 

‘W-l^  {^)''  (ii-o 

where  and  + Ao  ^ are  the  grasing  angles  at  the  source  of  the  pair  of 

rays,  is  the  nean  angle  at  the  receiver  and  Ar  is  the  horizontal 

••PMJation  of  the  rays  at  the  receiver. 

if  the  receiver  is  bracketted  by  two  rays  which  are  alaost  horizontal 
(i.e.  near  a ray  vertex  where  sin  approaches  zero)  the  following 

expression  is  used: 


coso^ 

1 /l  = . - 

r o rA  z 


Aa  a 


(II.2) 


where  A z is  the  vertical  ray  separation. 

In  the  vicinity  of  a caustic  the  energy  at  the  caustic  is  smoothed 
between  the  two  rays,  which  bound  both  the  receiver  and  the  oanatic. 

II. 2 Range  ealculations 

Expressions  will  be  given  for  calculating  the  horizontal  range  (r) 
travelled  by  a ray,  which  passes  downwards  through  a layer  of  the  ocean 
with  valties  of  depth  and  sound  speed  given  by  (^  , c^  ) and  (^ , ), 

where  Sj  is  greater  than  S| , and  is  greater  then  gS|  . Using  appropriate 

changes  of  variable  the  ease  where  greater  than  c^  and  the  ease  where 

the  ray  turns  over  within  a layer  (o,  than  beooaes  equal  to  o.)  can  be 

handled  by  the  expressions  given  b alow. 

The  range  ia  given  by  (e.g.  ref. 12  page  54) t 


uj  (1  - i*(s)/c?^)^ 


da 


(II.3) 


where  is  the  sound  apaad  at  the  ray  vertex. 


VSRL-0035-TM 


- 10  - 


Replacing  the  variable  r can  be  rewritten: 


dz 


where  Lt  - 


z.  and  z 


g - Zj 


(II.4) 


11.2.1  Linear  se^ients 

For  the  case  of  linear  aegnents  in  the  profile,  the  profile  in  the 
layer  become: 


C^C*  - 1 = 1 - Bz  (11.5) 

Substituting  11.5  in  11.4  the  range  in  the  layer  becomes: 

r = 24j  - (a-B)M  (II. 6) 


[0'  - -] 

Az 

Using  Oausslan  quadrature  (or  Oauaaian  integration)  the  intregral  oan  be 
evaluated  as: 

i - 1 

The  ooeffioiente  W.  ar^  given  for  a paztioular  value  of  n (e.g.  ref. 13) 
and  the  (z^)  values  am  given  by:  z^  ■ e ^ where  is  the  i th  zero 

of  the  Legendre  polynomial 

The  programme  has  included  aa  data,  values  of  and  z^  for  a equal  to 
3*  An  early  investigatioii  eittried  out  with  the  prngrzmme  ehowed  that  this 
order  gave  adequate  aoouzaoy  for  the  range  oaloulation. 


where  A » - 1 

2 2 

and  B 

C <3 

''  1 ''2 


11.2.2  Spline  fit 


Set  y 


Equation  II. 4 become:  r 


I 


eRA2mG  ANGLE  (•) 


I 
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